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Abstract 

We obtain an hybrid expression for the heat-kernel, and from that 
the density of the free energy, for a minimally coupled scalar field in 
a Schwarzschild geometry at finite temperature. This gives us the 
zero-point energy density as a function of the distance from the mas- 
sive object generating the gravitational field. The contribution to the 
zero-point energy due to the curvature is extracted too, in this way ar- 
riving at a renormalised expression for the energy density (the Casimir 
energy density). We use this to find an expression for other physical 
quantities: internal energy, pressure and entropy. It turns out that 
the disturbance of the surrounding vacuum generates entropy. For (3 
small the entropy is positive for r > 2M. We also find that the internal 
energy can be negative outside the horizon pointing to the existence 
of bound states. The total internal energy inside the horizon turns 
out to be finite but complex, the imaginary part to be interpreted as 
responsible for particle creation. 

1 Introduction 

One of the most important quantum physical quantities is probably the zero- 
point energy. It is needed, for instance, when we want to study back-reaction, 
i.e., the influence the matter fields moving in a curved back-ground assert 



on the back-ground geometry itself. This would be done by solving the Ein- 
stein equations with the expectation value of the energy-momentum tensor 
as source: 



where we have chosen units such that n — H — c — 1. It is also important for 
the study of renormalisation properties of a quantum field theory in a curved 
space-time. 

In order to find this quantity, we evaluate the heat-kernel corresponding 
to the equations of motion for a scalar field. This calculation is carried 
out at finite temperature. The integral of this heat-kernel with respect to 
some fictitious fifth coordinate, a, will give the density of Helmholtz' free 
energy (its integral over the entire space-time will give the zeta-function, 
which is essentially just Helmholtz' free energy). The resulting expression is 
regularised and renormalised in the subsequent sections. 
We do this by using the spherical symmetry of the spacetime in order to 
collect all the unknown bits of the heat kernel into a function g n i(r,r';a) 
depending only upon the radial coordinates and a. A recursive relation for 
an asympototic expansion of this unknow function can be found and solved, 
thereby allowing us to find the heat kernel. 

From the freen energy one can derive expressions for the entropy and the 
pressure by using standard thermodynamic relations. We find these and 
comment on their meaning. We also show that, to lowest order in the mass 
generating the Schwarzschild geometry, the free energy is that of an infinite 
family of particles moving in one dimension with an r-dependent mass. 



We consider a minimally coupled scalar field moving in a Schwarzschild 
back-ground, hence the action is 




2 Set-Up 




(i) 



where the metric is given by the standard expression 
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where M is the mass of the (classical) massive object generating the Schwarzschild 
geometry, a black hole or a star, say. We will put the mass, /i, of the scalar 
field equal to zero. We can later reinsert it if we find it desirable. 
The d'Alembertian becomes 

□ - -— - -—r 2 h— - — 



h dt 2 r 2 dr dr r 2 



where h(r) = (1 — 2M/r) and L 2 is the square of the angular momentum 
operator. A finite temperature can be included by complexifying the time 
coordinate t, t i— > r. As describes a Bose field, it becomes periodic in r, 
<f)(r + j3) = <P{t). This implies that the time direction is not only complexified 
but also compactified - by appropriate scaling we can take r to lie on the 
unit circle S 1 (see e.g. Ramond [|IJ or Itzykson and Zuber, 0, for further 
details on this). The action, furthermore, becomes "euclideanised" . 
Now, the partition function 



is simply a Gaussian and hence the functional integral can be carried out, 
the result being 

Z = (detD)- 1/2 . (2) 

See for instance [|]]. The quantity we are particularly interested in, is Helmholtz' 
free energy, F, which is defined as 

F = ~lnZ= JLlndetD. (3) 

The internal energy and pressure, which appear in the energy momentum 
tensor, are related to F by the usual thermodynamic relations, as is the 
entropy. 



3 Functional Determinants 

The major problem is apparently the calculation of det □. For consistency, 
and in order to establish notation, I will give a short introduction to the topic 
here. Descriptions can be found in e.g. Hawking or Ramond [|TJ. 
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A priori, the determinant of an operator A must be given by the product of 
its eigenvalues A 

det A = Yl A. (4) 

Now, obviously this is not an easy thing to calculate. This is where the zeta 
function, ( A (s), an d the heat-kernel, G A (x, x'; a), comes into play. Define 

caoo = ;>>-*, (5) 

then 

detA = e-&<°>. (6) 
The heat-kernel is defined through the differential equation 

d 

A x G A (x, x; a) = -—G A (x, x ; a), (7) 

subject to the boundary condition 

lim G A (x, x'; a) = 8(x — x'). (8) 

The reason for the terminology is transparent: when A is J^-, the Laplacean 
on S 1 , then Ca(s) oc C( s )) where ((s) = Yl^=i n ~ S i * s ^ e Riemann zeta- 
function, and when A = V 2 and a is the temperature, then the heat kernel 
satisfies the usual heat equation. 
Denoting the eigenfunctions of A by ip\ we have 

G A (x, x'; cr) = J2 MxWx(x')e- x ° , (9) 

A 

and one easily proves the important relationship 

Ca(s) = ^ J™ daa s - 1 J d 4 xG A (x,x;a). (10) 

Notice that the integral is only along the diagonal x = x'. We don't have 
to know the eigenfunctions in order to solve the heat equation, in a d- 
dimensional Euclidean space, for instance, one can show [|TJ 

G A (x,x';a) = (Anay^e-^ 1 , 
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in Cartesian coordinates. We will later use a mixture of techniques to arrive 
at a workable expression for the heat kernel. The angular and thermal part 
will be handled by a mode sum, whereas we will have to make do with an 
asymptotic expransion for the radial part in order to find the Casimir energy 
density. 

We can arrive at an important interpretation of the value of the heat kernel 
along the diagonal by noting 

is the density of Helmholtz' free energy 

F = J f(x)y/\g]d 4 x, (12) 
and is thus exactly the quantity we are interested in calculating. 



daa s Ga(x,x;ct) 



4 Finding the Heat-Kernel 

We now return to our particular problem, the minimally coupled scalar field 
in a Schwarzschild back-ground, i.e. the calculation of the determinant of the 
d'Alembertian. For this particular operator we know that the eigenfunctions 
can be written 

Mr,r,n) = e-^ T Y lm (n)g x (r), (13) 
where Q denotes the angles, and where cu n is the Matsubara frequency 

Un = — n = 0,±l,±2,.... (14) 

We have no intention of finding the eigenvalues A, nor of finding g\(r), instead 
we will rewrite the d'Alembertian as 



to. 



2 



1 d 2 d 1(1 + 1) 



D = -f- ? ^-V^ < 15) 

And the heat-kernel will be written as 

G(x,x';a) = ^^ Bl (r,r';(7)y lm (n)l^(n') e -- x -'^- <w -^, (16) 

nlm 
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where we have defined 

X m{r) = + y (17) 

The reasoning behind this formula is as follows. The eigenfunctions can be 
written according to fll3|) as a product of a spherical harmonic, a plane wave 
invloving the Matsubara frequency and the complexified time, and finally 
some unknown radial function. Clearly, we would expect the heat kernel, 
being a sum of products of eigenfunctions, to have a similar form, but since 
we cannot find the radial eigenvalues A we cannot simply use ([1"3| ) to find 
G. On the other hand, it must be possible to expand it on the functions 
Yi m (Q)Y l l l (Q') and e ~ ia, ™(' r ~ T '), hence the expression fllEp. 
A mass, /i, of the scalar field can be reinserted simply by adding jj 2 to X n i(r). 
The heat equation determines the unknown functions g n j. We will Taylor 
expand these in a, writing (asymptotic expansion) 1 

1 , 9 00 

g nl (r, r'\ a) = —j==e ~ \ a k (r, r')a k , (18) 

where the first term is the form g n \ would have in a flat space-time. The 
boundary condition implies a = 1. This asymptotic expansion is similar 
to the standard Schwinger-DeWitt expansion |5|, |6j. The difference here 
lies in the fact that we only make an asymptotic expansion of the radial 
part of the heat kernel, whereas we use a mode sum for the remaining part. 
Our expansion is thus a hybrid, combining the asymptotic Schwinger-DeWitt 
expansion with the exact mode sum. 

Inserting this expansion into the heat equation yields a recursion relation 
for the afc's. As we only need to know the value of the heat-kernel along 
the diagonal x = x', we will put r = r' to arrive at the following recursion 
relation by straightforward computation 

M\ 

l + k^ ctfe+i = -L x a k - £ 2 a fc _i, (19) 

r / 

where we have defined the operators 

2M\ d 2 (2 8M\ d 
Li = 1 ^7 + 



dr 2 \r r 2 I dr 



1 For simplicity we will follow the common abuse of terminology and refer to this ex- 
pression as "the heat kernel" although it strictly speaking is only an asymptotic formula 
valid for a not too large. 



6 



= h ^ + l (2h - 1 4r' < 20) 

u = h\' nl -^ + 2h\' nir - 1 + h'x„ l + h\'; d . (21) 

Putting 

H = L 2 a = 2h\' nl r- 1 + h'X' nl + h\'^ (22) 
the first few functions a& becomes 

a = 1, (23) 

ai = 0, (24) 

a 2 = (25) 

- = (») 

(~1) 4 2 4 T 1 1 H , (~1) 3 2% 1 

5 — "Tl I — Iri r If? r 1 r 1 r ~ TZ. ', — J-'l'p. V 



11-h 9-h 7-h 5-h 11-h 9-h 5-h 
-1) 3 2 3 „ 1 , H 



, L 2 -L x -. (28) 

11-h 7-h 5-h K ' 

There are singularities in these at r = 2M, as one would expect. Noting that 
both H and L 2 are linear in lu^, we see that can contain all even powers 
of the Matsubara frequencies up to and including w„ fc ^ 2 ' . 
The method of summation of modes when calculating zeta functions have 
been applied in a number of other spacetimes by various authors, [10|, and 
a number of explicit results exist for some rather simple cases. 



5 The Free Energy Density 

We have now found an asymptotic expansion for the heat-kernel and are 
ready to integrate it. The a integration is trivial (it usually is), and will just 
give a sum of T-functions. Explicitly 

1 /"OO 1 °° 



7 



JL Eai(r) £(l±±J) A | i — . (29) 

nZm fc=0 



One should note that even though the asymptotic series is only valid for o 
small, this integral still makes sense since the factor \ n \ acts as a mass term 
regularisation. If we hadn't used our hybrid method we would have had to 
assume a mass for the scalar field in order to obtain a meaningfull integral. 
We can still do this, if we wish, by simply adding fi 2 to A n ;. The addition of 
such a finite mass, however, will not modify the convergence properties of /, 
hence the succesive regularisations to be performed below are still needed. 
This is so because only the sum over n need to be regularised, and for fixed 
I, X n i acts like a finite mass irrespective of whether /i 2 7^ or not. 
For k 7^ 0, differentiating this with respect to s at s = simply amounts 
to removing the T-function in the denominator and putting s = in the 
remaining terms. The k = term is singular and has to be treated separately. 
We have 

ft \ — IV |2 1 ~~ 2' \"l~ s 

fk=0{X; S) = \yim\ -7= T , j Kl , 



nlm 

so we have a divergent sum (for s = 0) in the numerator, namely V z 2 + Aix 2 n 2 
where z 2 = K ; +i)fc(r)/? _ ^y e can carr y ou t the sum for an arbitrary s, see e.g. 
Ramond Iffl: 



\ 00 1 
— J2 (^ 2 + 47rV)^ 



r ( s ) ^ v > r ( s ) 

v « n=— 00 v « 



1.2. , 2C(2a+l) 

•2 



(4vr 



2 ^x 2 ^ (-If , L , 

— r-> — i ^ — r r(s + i/--)C(2s + 2i/-l 

Differentiating and putting s = we obtain 

1 2 o 00 2i^ (1^ 1 

z + ^(7 - 2 ln27r)(47r - -) + — £ £--LJ r(i/ - -)C(2«/ - 1), 

which is finite. This expression is to be multiplied by /J" 1 /^ 1 / 2 . 
The free energy density f(x) is then 

f( \ 1 d f 

f{x) = ~W^s 

s=0 



8 



i in 



1 



-2^ 



1(1 + 1) 



2 v^^(7-21n2vr) 



4tt Z(Z + 
Vh(3 rrr 2 J 



+ 



1 °° 1 

"og E l^| 2 ( 471 ") _1/2 E « fc (r)A 7l z(r)^ fc r(A; - -). 



(30) 



On the face of it, this function depends on all the coordinates, but looking 
closer we see that there is no explicit time-dependence, and using the sum 
rule 

21 4- 1 

m 

where u, u' are unit vectors given by the solid angles Q, Q', we see that H 



47T 



21 + 1 



E 



'-lf(2l-2v)\ 



^ ^y\(l-y)\{l-2v)\- 



So the angular dependency also disappears. In accordance with what we 
would expect, the energy density depends only on the radial coordinate. We 
will write the density as a mode sum 



(31) 



where 
fi(r) 



vr -l^-2_L( 47r )-3/2( 2/ + 



1(1 + 1)(3 2 



2vr 3 / 2 



(7 - 21n27r) 



/ 4tt /(/ + l)V^/9 2 

\7W 2 
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2 eW^)V->- 



4vr 2 r 2 7 v 2' 



x 



1 oo oo ^ 

^(4vr)- 3 / 2 (2/ + l)2-' E r (*4 )fli{rWr)VW 

n=— oo fc=2 J 

Ijg (-l)»(2I-2i/)l 

A.^-v)!(Z-2z,)!- 1 J 

The major problem here is that the coefficients eifc(r) depend on n through 
the Matsubara frequencies. 

This expression needs regularisation. For instance, when I = we encounter 
a singularity. Remembering that a k can include all even powers of n up to 
and including r^ k l 2 ^ we write 

m 

a k {r) = Y J bf{r)n 2 \ 
t=o 

where b k s ' is independent of n. Defining A = Ait 2 /hj3 2 we then have 

oo oo oo .. [^/2] oo 

n=-oo fc=2 fc=2 t=0 n=l 

[fc/2 



oo 1 [K/-4J 

2^^- ^)A^ fc £ bf (r)((2k - 2t - 1) 



fc=2 t=0 



which is divergent whenever 2k — 2s — 1 = 1, as s < this can only happen 
when fc = 2, s = 1. 

6 The Minkowski Space Contribution and the 
Casimir Energy Density 

To get a physical quantity we must subtract the Minkowski space contri- 
bution || ||, because the physical energy must be such that it vanishes in 
flat spacetime. Doing this we thereby get the "Casimir energy density" . The 
relationship between this Casimir energy density, (^-functions and cut-off reg- 
ularisations has been studied in 11. This subtracting off the flat spacetime 
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contribution is sometimes enough to cancel divergences, because any man- 
ifold looks locally like Minkowski spacetime and hence will have the same 
leading divergence as in flat spacetime. It will turn out, however, that the 
Casimir energy thus obtained in this case is not finite and still needs some 
regularisation. The subtraction off the flat contribution should therefore not 
be seen as much as a regularisation/renormalisation as simply a normalisa- 
tion, 

The heat kernel for the d'Alembertian is already known in flat spacetime 
using Cartesian coordinates, but this will not be useful to us here, since we 
need to act on the asymptotic expansion Hence what we need, is to 

find a similar expression for the flat spacetime heat kernel using a mixture of 
mode sum (angular coordinates and Matsubara frequencies) and asymptotic 
expansion (radial coordinate). We obtain this contribution by letting M — > 
in the metric (and hence the d'Alembertian etc.).Q Denoting the resulting 
coefficients by dk(r) we get the following much simpler recursion relation 

1 ? „ 1 ~ _ 
k + l k+l 
1 ( d 2 2 d \ 21(1 + 1) ( d 2\ 



k + l \dr 2 r dr J (k + l)r 3 \dr r 
The solutions go like 

d k (r) = a k (l)r- 2k , (34) 

there is no dependency upon the frequency uj n . The coefficients a k satisfy 
(for k ^ 0) 

2fc(2fc-l) 4kl(l + 1) 

a k+1 = — at+i ak + nrn - "*- 1 ' (35) 

and ao = 1, «i = 0. Written down explicitly, the first few of the d k areQ 

d = 1 = a (l), 



(ii = = ai(/)r" 



-2 



2 This is not the same as saying our result for the heat kernel is smooth in M. We merely 
use the M — > as a short hand for saying that in that limit the metric tensor reduces to 
the Minkowski spacetime one - with a compact time due to the finite temperature. The 
coefficients a k will not a priori be smooth in this limit, however. Hence, new recursion 
relations have to be found. 

If the mass, [i, of the quanta of the scalar field is non-zero, then these would be com- 
pletely different. They would still be independent of the Matsubara frequencies, though. 
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a 2 = 2l(l + l)r- A = a 2 (l)r- A , 

a 3 = -8/(/ + l)r- 6 = a 3 (Z)r- 6 , 

5 4 = 12/(/ + l)(5-/(/ + l))r- 8 = a 4 (0^ _8 > 

a 5 = -672Z(Z + 1) - -^1(1 + 1)) r- 10 = a 5 (Z)r- 10 . 

Denoting the corresponding free energy density by f(°\ the Casimir density 
is defined to be 



/ Cas (r) = /(r)-/(°)(r) = ^/ z Cas (r). 



(36) 



The very simple form the functions take on actually allow us to carry out 
the summation over the Matsubara frequencies, using essentially the same 
formulas as in the M^O case. Defining 



[1/2] 



Cl 



1 fA^-V2,o, , (-l)"(2f-2i/)! 

¥ 4n) (2/ + 1)2 hw^W^y. 



(37) 



the resulting free energy becomes 



n 



Cas 



2-v/7TQ 



7 -21n2 V/i _ 1/2 _, 2 _^ + 4^ )+ 



2tt 5 / 2 



v^/3 2 /9 2 



(3^ 



X^-V _ ^) ( 1 Y{u - ±) C (2,/ - 1) + Y: T(k - h x 



v=2 



Z/!(47T 



fc=2 



(I) 



Y: a fe (r)Ajr fc -« fc (/)r- 2 ^ 2fc - 1 (2 7 r) 1 - 2fc 

1 = — oo 

o 00 ~2u /_i \i/ 1 

V(k - \) ^ i/! (2tt) 2 ^ v 2 ;sv 



l-2k 



where 



^ 2 _ Z(Z + l)h(r)p 



2 = Z(Z + 1^_ 2 



+ 2C(2A;- 1) + 
(38) 

(39) 



and where we have used zh' 1 ^ 2 — 5 = 0. We should note that con- 
tains a countable infinity of divergences, this time coming from the terms 
z l ~ 2k a k (l) ~ 1(1 + i) 1 -2fe+[fe/2]_ We should furthermore notice that all the 
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terms z 2v hr x l 2 — z 2v are finite as r — * 2M, since zhr 1 ! 2 — > in this limit. 
The only divergent parts as r — > 2M are the afc(r)A2 _fc -terms. We can find 
an improved expression for these by writing once more 



t=o 

we then have to perform sums of the form 

oo 

Us,t;a;z)= ]T n 2t (z 2 + an 2 )" 5 , (40) 

n=— oo 

with a = 2tt. For £ = we know the result; it is just the formula we have been 
using a number of times by now. Denote this sum by £i(s; a; z). Explicitly 
it is 



fi(s;a;z) = (z 2 + an 2 ) 



2\-s 



9 ( — ■\\ v 7 i v 

r(s) ^ z/!a s+iy 



For t a positive integer (which is the only case we need to worry about) we 
notice the following relationship 

6(s, t; a; z) = L^L-y — _ t . a . z y ( 41 ) 
Inserting the expression of £i we finally arrive at 

v=0 



6( ., i; a; 2) = _Ei£L E( _^£l£|l±| c(2s _ 2( + 2 „ )a -.-«-, (42) 



We have to evaluate this at s = k— |, a = 47r 2 and multiply it by ({3 2 h) k 2 2y / 7TQ 
in order to get the contribution to the free energy, which then reads 



/, Cas (r) = 2V5Fq 



7 ' 21n27r f^- 1/2 - z 2 - — + —) + 
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2 °° 



!/=2 

oo / [k/2] , _ co 

E r ( k - ^r^l 2 E E(/^) M (-ir 

fc=2 "' y t=0 ^ 2 / ^=0 

F(fc - I - 1) ( 27r )~ 2fe ~ 2 " +1 ~ 2 *C(2^ + 2A; - 2* - 1) - 

l-2fc 

a fc (/)r- 2fc /3 2fe - 1 (27r) 1 - 2fe <! ( — ) + 2((2A; - 1) + 




F(*3TyE^ F r( 4 + ,- 5 )c(2* + 2,-i) 



We note the presence of a singularity when fc = 2, t = 1. It is the same kind 
of singularity we encountered in the / = case. The regularisation of such 
singularities is the subject of the next section. 

It is an easy matter to derive a recursion relation for the new coefficients b^; 
inserting = ^2 t b^n 2t into the recursion relations for the yields 

(i + * + e = -Mi? - if e, - W, (44) 

where we have written L 2 = + n 2 L^\ Thus 

4»» = - 2 >!±ll(4 + h '- h l' 

r 6 \ dr r 

4» = -^CvA+fli-I-^ + fc- 



/i/3 2 \ dr r h 

The first few of these have been written out explicitly in the appendix. 
Remembering that vanishes whenever t < or t > [|] , we can calculate 
these coefficients systematically. We furthermore notice that is (3~ 2t times 
some function which only depends upon the metric and its derivatives (i.e. 
on h, h', h", ...), we can thus write 

b®(r,f})=r*<$\r), (45) 

in this way separating the coefficients into a purely thermal and a purely 
geometric part. 
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Now, from Helmholtz' free energy density / Cas we can calculate various in- 
teresting physical quantities, namely the renormalised expressions for the 
(modes of) the pressure, pi, internal energy density, ui, and finally the en- 
tropy density si by 



dV J g Arrr 2 \ dr 
ui = pff» + p 8h (47) 



Cas 



However, if one attempts to calculate, say, the entropy density, one will find 
that it is divergent; thus there is still some regularisations to be done. 



7 The Final Regularisations 



As we saw, / Cas , had a singularity from 1 = 0, stemming from two kinds of 
singularities in / and which did not cancel each other. We furthermore 
noticed a similar singularity in the k = 2,t = 1 term. We will now return to 
this problem once more. 

The singularities in the curved and flat space-time free energies came from 
two different sources, one was C(l); i- e - the pole of the Riemann ^-function, 
whereas the flat space-time contribution had a countable infinity of poles 
(1(1 + l))~ v - We thus need two different ways of dealing with the problems. 
It is well-known that singularities in the ^-function regularization are to be 
removed by taking the appropriate principal parts (see e.g. ||), which for 
a meromorphic function simply amounts to extracting the finite part near a 
pole. Hence 

-3/2 



Mr) 



-V7T 



4vr 7 



1=0 



(r)C(3) + 7 & 



1=0 



k=3 



where we have written 



4tt2 



t=0 



(t) 

k 



1=0 



(r)((2k-2t-l) 



[k/2] 

a k(r) = ^bf(r)n 
t=o 



21 



(49) 



(50) 
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as before. In obtaining this result we have used |fj. [7J 

^ oo 

«*) = — 7 + 7 + E(- 1 ) n2 r( s - 1 )"- 

S — 1 TV. 



n=l 



Now 



-2H 



a 2 



5-h 

(2T- 1 -h-\h'f + h") (for 1 = 0), 



^n 2 n 2 



so fef 



i=0 



(5 - h)h(5 2 
leading to 



/o(r) = Y ^Vh 1 (5-h)-\2r- 1 -h-\h') 2 + h") + 



i ^ i / /T^\ 2fe-1 [fe/2] 



where we have inserted c = — ^g(47r)~ 3 / 2 . We notice the appearance of a 
temperature independent term, this will not give any contribution to the en- 
tropy then. 

The remaining singularity in fi for / > 1 also comes from the pole in Rie- 
mann's zeta-function and again appears only in the k = 2 contribution. Here 
gets multiplied by 



00 ? 2v o 

A^c^tf' 2 ^(-l)--r(z/ + -)C(3 + 2i/)(4tt 



which is non-singular, but b 2 gets multiplied by 

4^FQ^ 3 / 2 (C(l)r(^)(4vr 2 )- 3 / 2 + ...), 

with the non-singular terms left out. In the spirit of ^-function regularization 
then, we must interpret this £(1) as 7, and the result is then well behaved. 

4 The n = contribution will always be proportional to 1(1 + 1), as this is the only 
w n -independent term on the right hand side of the recursion relations. For k = 2 we find 
explicitly = Ar~H{l + 1)(5 - h)- l {ti - h/r). 
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For the Minkowski space contribution we have to return to the initial prob- 
lem, the solution of the heat-equation, and we have to redo the calculation 
with 1 = 0. Hence 



/o [r) ~ 2(3 ds 



s=0 




(52) 

where solves Dg^ = w*g$ - £g§o- So 9n0 is the Minkowski space 
analogue of the function g n i we introduced in order to solve the heat-equation 
in a Schwarzschild geometry. The solution is easily seen to be 

g n Q = e ^~ + terms vanishing at r = r'. (53) 

This leads to (standard calculation) 

We notice that this contribution is independent of r and that it vanishes as 
P — > oo, it is thus a pure effect of the finite temperature. This also implies 
that there is no Minkowski spacetime contribution, for I = 0, to the pressure, 
but only to the entropy and internal energy. It is interesting to point out 
that in fact the entire Minkowski space contribution is /3-dependent, also for 
1^0. 

The I = contribution to the Casimir energy density can now be written 
down: 




The temperature independent part of this is negative for r not too far away 
from the Scwharzschild radius, and in fact looks pretty much like the usual 
effective potential (i.e., Coulomb plus angular momentum part) for the hy- 
drogen atom, thereby suggesting the existence of bound states. We will later 
see more convincing arguments for this. 
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From f^ s we can evaluate the / = contribution to the entropy density, 
which turns out to be 



1 00 1 / fha\ 2k ~ 1[k/2] 



fc=3 



i=0 



91tt/3 

(56) 



For I ^ we still had a singularity, using the above prescription for its 
regularization we arrive at (including the / = contribution) 



fp as (r) = 2^Q 



7-21U2. / _ 1/2 _ _ 
2tt 5 / 2 V 



4tt 2 4tt 2 ' 



E 



ir - 1 {zV-~z^)Y{v- W-l) + 



/? 2 

1 

2 



i/=0 



v 2^ 



^ (4^^/V-v 2 g g^r(§ + , )C (2, + 3)+ 

4 1) / i 3 / 2 (2vr)- 5 
a 2 (/)r- 4 (27r)- 3 



_|_ r( , + 3 ) _L c(2l/ + 1) 



2,r\ 3 4 



^=0 



1 / [fe/2] m-i 
yr(k- V^ 1 2V Ufc . 2 



;-i)^ 



r(fc - 1 + 1/) 
r(fc - i-t) 



2 7 C(2A: - 1 + 2i/ - 2f) - a k {l)r' 2k {2n) l ' 2k 



z/!(27r) 2fc+2iy+2 ' _1 

l-2fe 



X 



r^^r £ L T7^r(A; + 1/ - -)C(2* + 2i/ - 1) 



v^7 



2?' 



+ 



(2r- 1 -/i- 1 (/i , ) 2 + ^") + 



16tt 3 / 2 (5-/i) v "' v "' ' " ' ' 91tt/? 2 

This leads to the following expression for the modes of the entropy density 



(57) 



7T 
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2 y (- i n^- 1 ) {z 2, h -i/2 _ ~2 V)v{v _ i )c(2 _ x) 

v^F^ i/!(2tt) 21 '- 1 V ' V 2 ;sv ; 

i^/j* g ti£g^) r (, + |) { (2, + 3)+ 

v ff» ( 9hk -i/2 v r ( fe -D h w v (-ir^+2t+2fc-i) 

fc=3 y *=o v 2 / ^=0 v I 

T{k-\ + v 



\ '- ((2v + 2k - 2t - 1) - 



Y{k-\-t) 

2*^1-2* 

r(fc - |) ^ ^!(27r) 



2a fc (/)r~ 2fc (2 7 r) 1 - 2fc - 1 lN V ( ./f/ r^ + f - -) x 



(2k + 2u- l)C{2k + 2u- 1))] - (58) 

One should note that is negative in some regions, this need not indicate 
that the entropy as such is negative, but only that one can not really localise 
the entropy: the correct physical quantity is S — J2i I Sir 2 dr and not Si(r), 
and this could very well be positive still. Furthermore, the entropy density 
is not uniquely defined as one could add a total derivative and still get the 
same over all entropy. This arbitrariness in the choice of si will of course be 
fixed by defining a zero for the entropy proper. Furthermore, the calculation 
is only valid outside the horizon, r > 2M, so in principle one could have a 
negative entropy in this part of the system (here the universe) , provided that 
a suitable positive entropy is present inside the horizon r < 2M. It does 
seem rather strange, however, to insist on this interpretation, as r > 2M 
is the only observable part of the universe, and we would certainly expect 
physical quantities to "behave properly" in this region. 
The entropy density we have found here is partly geometric in nature (in- 
duced by the curvature) and partly thermal (coming from (3 7^ 00). The "ge- 
ometric entropy" of black holes has been studied by Moretti in |ll] , where he 
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is also using ^-function techniques. He shows that the Bekenstein-Hawking 
entropy is purely geometrical. We refer to his paper for further details. 
The /3-independent part of / ; Cas is seen to be simply (remember, q ~ b 



tf"L t em , = ^c^h^r 1 -, + y J^_ (2r- 1 - h-\h'f + h") 



2W/ 2 1 ] ^r + M) 2^ ul{l _ u y {l _ 2u)l + 



-i^v^ 7 (5 - /^)- 1 (2r- 1 - ft" 1 ^') 2 + h"), (59) 

1d7T ' 

(the last term coming from I = 0) which is divergent on the horizon. By 
using 

~' 2/ + 1^1 (-1)"(2Z - 2u)\ 



Ely |2 _ 1 ^ 2/ + 1 ^ 



47r ^ 2' ^ v\(l - u)Ul - 2u)V 

Ira 1=0 u=0 v ' y ' 

which follow from the formula for Pi(l) used before, we can actually carry 
out the sum over / to obtain 

/ Cas | t = ^-r- 2 (r-2M)-^ 2 (2r+M)- l + ^r^-hrHh'Y+h"). 

J Inotcmp 167r 5/2 V J V t , t ^3/2(5 _ ^ V V ) ~r , 

(60) 

It is also interesting to note that outside the Schwarzschild radius this is 
positive, whereas as inside it is negative imaginary. The integral, moreover, 
over r is finite (for / 7^ 0) in any case, even though the function is divergent 
for r = 2M, and we get the following result 



7 



temp lQ n 5/2 



71 / 2 ,/2\ 

— 1>\ ^rrArtanh —= 



VTOM V 5M \y/EJ 



I > 1, (61) 



the first term is the r > 2M contribution, whereas the second comes from 
< r < 2M. This result represents, then the energy coming from the 
removal of the point r = from Minkowski space and placing a mass M in 
that singularity. It is, in other words, a proper Casimir energy similarly to 
the one one obtains in the original case of two plates in flat spacetime. 
From (|60|) we can also find a temperature independent pressure by simply 
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taking the derivative with respect to r. Doing this we find 
M 2 7 8M 2 + 19Mr-Ur 2 



P\ 



no temp 2 r 4 /l 3 / 2 (2r + M) 2 

7 10M 4 + 16M 3 r - 32M 2 r 2 - 2M 3 r 2 + 12Mr 3 - HM 2 r 3 + lOMr 4 



64tt 5 / 2 



r 5/ l 3/2( 2r + M) 2 



This local pressure is negative outside the horizon, i.e., for r > 2M, and 
negative imaginary for r < 2M. 

For small but non-vanishing f3 we get (writing q = /9 _1 q,z 2 = /5 2 Z 2 ,z 2 = 
(3 2 Z 2 ) 



ft 



Cas 



2-v/vtq 



7 - 21n2vr 

2tt 5 /2 



P(Z 2 h-V 2 -Z 2 - (3-Htx 2 ^ 1 ' 2 - 1)) + 



V 7 / 2 /^ 2 /*" 1 / 2 _ Z 2 )r(3/2)C(3) + 7 ^4 0) / i 3 / 2 /? 2 r(3/2)C(3) + 

8 27T 

(2vr)- 5 / i 3 / 2 4 1) Q 7 " (27r)- 2 Z 2 /? 2 r(5/2)C(3)) - 



a 2 (/)r- 4 (27r)- 3 ( (27r) 3 Z" 3 /?- 4 + ^=/? 2 r(3/2)C(3) 



oo 1 

^r(A;-i)a fc (/)r^Z 1 - 2 ^- 1 



fc=3 



+ 



v/ '" : .(2r~ 1 -h~\ti) 2 + ti')+ ' 



16tt 3 / 2 (5 - /i) 



(63) 



Note that in the summation over fc contribution, the summation does not 
extend to a summation over powers of /3; is multiplied by a function 
of r,l, A,(r) := J2T=3 F ( k ~ l/2)a fe (0 r ~ 2fc 'Z x ~ 2k . The free energy resulting 
from this expression turns out to be everywhere negative for r > 2M and 
imaginary for r < 2M. 
The entropy density we derive from this is 



^Cas 



= 2a/7TQ 



'- ^- 2 t 2 \ z 2 h^ 2 - Z 2 ) - 6r 2 7 -21n 27 r (/rl/2 _ 
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-^Cm\Z 2 h- l l 2 - Z 2 ) + i 7 r- 1 /2 c (o)/ l 3/2 C(3)/3 3 _ 
8 2 

^(27r)- 7 Z 2 / i 3 / 2 c 2 (l)C(3)/3 3 + Aa 2 (l)r- 4 (Z- 3 (3- 3 - C(3)(2tt)- 3 /3 3 ) + A, (7 



- 2r 
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(64) 



This is positive outside the Schwarzschild radius. We see that the tempera- 
ture independent part of the entropy density is simply Ai(r) for I > 1 which 
is independent of the mass of the black hole and comes from the renormali- 
sation, and the I = contribution then has all the information about M. As 
well the free energy as the entropy and internal energy densities are divergent 
on the horizon r = 2M, although in a rather mild way. 
The internal energy density is found to be 



u 



Cas 



2v^Fq [-( 7 - 21n27r)7r- 5/2 / 3 3 (Z 2 /r 1/2 - Z 2 ) - 47r 1/2 ( 7 - 2\n27i)p~\h' 1/2 - 1) + 



PXWiZ 2 ^ 1 ' 2 - Z 2 ) + -^ c fc(3)/i 3/2 /5 4 + 



16tt 3 ^ swv ' 4v^F 

{2.)-W 2 ^\\ 1 p 2 + ^ 2 Z 2 a^)- 

a 2 (/)r" 4 (-3Z- 3 r 2 + ^C(3)/? 4 ) + 

'' /; ''' .(2r- 1 -h-\h') 2 + ti')p 2 - — + 0(f3 5 ), (65) 



167r 3 / 2 (5-/i) v v ' " 91tt 

which is seen to be singular at r = 2M as mentioned above. It is interesting to 
note that, contrary to the free energy and the entropy densities, the internal 
energy is integrable inside the horizon. In fact we get 

^ Cas | int = / rfrr 2 wp as (66) 
Jo 

= -12M 2 (/(Z + I))- 1 ' 2 ?- 2 + 7M1 ^ 5 (24tt 3 / 2 - 32 - 5Mtt 3 / 2 )/3 2 + 

128V57T 3 

M 3 /T 1 (3647v^F - /5 - 728 v^F In 2tt) - 



273tt 



/5- 1 (1280M 2 7r( 7 -21n27r) + 7 /3 3 (M-8)) 



128^ 

{m + ln5) 128^v/ 2(32 " 24?r3/2 + 5M7r3/2) ' (67) 

which is clearly complex but finite. The imaginary part of this internal 
energy is interpreted as giving a continous creation of particles. Notice, by 
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the way, that only one of the terms in Up as is /-dependent. The temperature 
independent part, moreover, is seen to be — pf^, a negative constant. In 
any case, for all values of (3 there is a region in which Up as < 0, the size of 
this r-interval, however, becomes rapidly smaller as f3 increases. The actual 
depth of the resulting potential "well", moreover, also rapidly decreases as 
j3 grows. But for (3 ~ .1 a very large range of r values exist for which the 
internal energy is negative, once more pointing towards the possible existence 
of bound states. This region only exist for I < 10, though. For higher values 
of angular momentum, no such region seems to exist, suggesting that only 
orbits with low angular momentum can be bound, which is what one would 
expect on physical grounds anyway. 

Inside the Schwarzschild radius both the real and imaginary parts of the 
internal energy density seem to be negative for all values of (3, with the 
energy density becoming almost entirely a large negative imaginary number 
near the horizon and a large negative real number near the origin. This 
apparently holds for all values of I. 



8 Interpretation of / 



Cas 



We can find a closed expression for the summations over powers of z, by 
noting with Ramond |1[ that 

9 i*±i roc 

V*(d-2)\\J 

Inserting this at an earlier stage (i.e. before the cr-integration) we arrive at 
the following closed formula 



, i>--)C(2z/-l) oc -f3~ l Ti-' i/2 I ^w 2 + y 2 +2 In (l - ^ dw, 

(68) 



u=0 

where 



2tt 2vrr 2 



Now, this is essentially the expression for the free energy of a scalar quantum 
field in d — 1 dimension (see e.g. [|lj) and with the energy of the individual 
quanta being given by us 2 = y 2 / (3 2 = 1(1 + l)h(r) / (27rr 2 ) . Our Casimir energy 
density is hence the regularised expression for the energy of an infinite family 
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(labelled by their angular momentum quantum number /) of such quanta. 
This was to be expected, since the spherical symmetry of the system im- 
plies that the problem is essentialy unidimensional (the radial coordinate). 
Moreover, the angular quantum number I appears as defining an r-dependent 
mass. 

Expanding to first order in M we then arrive at expressions of the form (valid 
only for 1^0) 



w 2 + - ^— + 21n|l-eV 

2nr 2 



dw 



~ M T^T fVv + / ft( i+ l)))-V2 x 







1 + 2 e ——= ( w 2 + ^ + i >f ) )dw + 0(M 2 ). 

I _ e -y/w*+i(!+i)p/ r * V 2vrr 2 J J 

(70) 

The first two terms are the contribution coming from a massless particle in a 
flat one dimensional space. While the second term shows how the pressence 
of a gravitational field (here M) modifies the energy. 



9 Discussion and Conclusion 

We have obtained asymptotic expressions for the heat kernel and thereby for 
the free energy density in a Schwarzschild geometry by use of the ^-function 
technique. We also took the M — > limit, i.e. the flat space limit, and we 
subtracted the two energies, to obtain what we called the Casimir energy 
density, this is the part of the zero point energy density due to the curvature 
(i.e. to the deviation from Minkowski spacetime) and is thus an intrinsically 
interesting quantity. As we would expect, these densities all turned out to 
depend only on the radial coordinate. 

The major unanswered question, however, concerns the boundary conditions. 
The asymptotic expansion of the heat kernel is independent of the chosen 
boundary conditions, i.e. on the particular vacuum state. It is known, how- 
ever, that the full renormalised energy-momentum tensor is sensitive to the 
particular vacuum state. Page, E3], has computed the energy momentum 



tensor in a Schwarzschild background for conformal coupling, and he got the 
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00 contribution to be 



T ° « (-9216M 10 - 21504M 9 r - 18688M 8 r 2 + 7168A/V - 
2560MV + 4M 2 r 8 - 4Mr 9 + r 10 )/(122880M W) 



Anderson et al, |13| have extended the original calculation by Page to arbi- 
trary coupling to curvature in the Hartle-Hawking vacuum state, and they 
then find a finite energy density at the horizon p(r = 2M) = (15£ — 
4)/(15360M 4 7r 2 ). They also find that p is positive for r > 2M only for 
jjr < £ < 1.2575, where £ is the non-minimal coupling. In our case we have 
£ = 0, and we found the Helmholtz free energy to be positive for r > 2M 
in the limit of vanishing temperature. We also found, however, the internal 
energy to be divergent at the Schwarzschild radius, and, moreover, to posses 
regions in which it was negative, suggesting the existence of bound states. 
Our computation seems to have more in common with Boulware vacuum. 
Candelas, [jl5|| , has computed the renormalised value of ((f) 2 ) in the Boulware, 



Unruh and Hartle-Hawking vacuum states, and found that for r — > 2M the 
Boulware vacuum expression diverges, whereas the others are finite. Fur- 
thermore, as r — > oo, (0 2 ) vanishes in the Boulware vacuum and not in the 
otehr two. Both of these features are reproduced by the present calculation, 
which therefore seems to be closer related to the Boulware vacuum than any 
of the other two known vacuum states. It is impossible to tell precisely, 
however, since the Schwinger-DeWitt expansion which we used partially in 
obtaining our result is insensitive to the details of the vacuum state. This 
does not imply, though, that the present calculation is meaningless from a 



physical point of view, since Anderson and coworkers, ||14|| , have shown that 
for (r — 2M) /M not too large, the Schwinger-DeWitt expansion expression 
for (T 00 ) is valid, and only for (r - 2M)/M > 1 (i.e., r > 3M) does the 
state-dependence become noticable. Hence, for (r — 2M)/M < 1 at least, the 
present calculation is valid. Moreover, since we use a hybrid technique mix- 
ing the exact mode sum with the asymptotic Schwinger-DeWitt expansion, 
one can expect the region of validity of this calculation to be larger than the 
pure Schwinger-Dewitt approach. It remains for future research to specify 
the precise range of validity, to see just how far away from r « 3M we can 
push the present computation. 

All the calculations were carried out at finite temperature. The quantities 
we have found are therefore also the effective actions for a free (minimally 
coupled) scalar field in these backgrounds. The Casimir density then shows 
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how the presence of curvature modifies the effective action for a flat space- 
time. 

We noticed, for instance, the presence of a constant (with respect to r), 
positive contribution to the entropy, which we can interpret as generation of 
radiation, somewhat related to the Hawking radiation but not quite the same, 
as this time the radiation clearly came from the distortion of the surrounding 
vacuum and could not be attributed to the internal (quantum) structure of 
the massive object generating the Schwarzschild geometry. This contribu- 
tion, moreover, depended upon the temperature (was, in fact, proportional 
to it) and can in this way be seen as a correction to the Hawking radiation, 
by letting the temperature be equal to the Hawking temperature T H . 
Especially interesting would be the coupling of the Casimir energy density to 
the curvature, i.e. plugging in the Casimir energy density into the right hand 
side of Einstein's field equations and assume M = M(r, t) is a sufficiently 
slowly varying function of r,t (if not, one will have to redo the entire calcu- 
lation of / Cas with a time and radial dependent mass). Will the singularity 
at r = 0, stemming from the point mass (M(r,t) = M8(r) independent oft), 
become "dressed" in this way? If so, then quantum effects could be responsi- 
ble for the removal of other singularities as well, most notably the singularity 
at the Big Bang (or the Big Crunch). Will this constitute a general "quantum 
removal of singularities" -mechanism? This back reaction would also give us 
a more precise picture of the internal structure of a black hole, as well as of 
its stability; we already know that black holes can evaporate due to Hawk- 
ing radiation, but the calculations put forward here suggests the existence 
of even more such effects, this time coming directly from the disturbance of 
the surrounding vacuum, from the "dressing" of the black hole so to speak. 
This means that at least part of the Bekenstein-Hawking entropy must be of 



a geometrical nature, consequently supporting Moretti's findings, [lT|. The 
fact that we get negative values for the candidate entropy density shows that 
this entropy cannot be localised in one particular region (such as near the 
Schwarzschild radius, say) but has to be considered a global object intimately 
related to the global topological properties of the spacetime manifold. 
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A The Explicit Form of the 

In this appendix we list, for the readers convenience, the explicit form of the 
first few For k = we have simply = 1 while for k — 1 we have 
bf^ = 0. All of these coefficients can be written as rational functions oir,l,(3 
and M. 
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A.0.1 The Case of k = 2 

Here we have two possible values of t, namely t — 0, 1 and the functions are 
readily found to be 

,(o) . ^ r-AM 



b^> = -21(1 + 1) 



2r 5 + Mr 4 



(1) 16vr 2 M 2 
u 2 — 



/3 2 r(r-2M) 2 (2r + M) 



We note the appearance of a zero at r = AM in 6 2 ^ and a pole of order two 



at r = 2M in 6 2 



(i) 



A.0.2 The Case of k = 3 

Again we have only two coefficients, corresponding once more to t = 0,1. 
The functional expressions, however, become somewhat more complicated. 

(o) „, . <s 96r 4 - 640Mr 3 + 156M 2 r 2 + 288M 3 r + 64M 4 



^ = Z(Z + 1) 



r 6 (2r + M) 3 (3r + M) 



, m 9 o 768r 4 + 512Mr 3 + 244M 2 r 2 - 576M 3 r - 128M 4 

3 /5 2 r 3 (r-2M) 3 (2r + M) 3 (3r + M) 

Here there are no zeroes but only a pole of order three at r = 2M. 



A.0.3 The Case of k = 4 

For k = 4 we get very complicated expressions for all three possible values 

of t, t = 0, 1,2. 

The numerator for t = is 

1(1 + 1){(8640 - 103680/(/ + l))r 9 + (995328 - 37152/(7 + l))r 8 M - 
(532800 + 53712/(/ + l))r 7 M 2 - (1407552 - 79016/(/ + l))r 6 M 3 + 
(71328 + 207464/(/ + l))r 5 M 4 + (1035336 + 182480/(/ + l))r 4 M 5 + 
(717648 + 85024/(/ + l))r 3 M 6 + (228192 + 22504/(/ + l))r 2 M 7 + 
(36096 + 3208/(/ + l))rM 8 + (2304 + 192/(/ + 1))M 9 } 

while the denominator is 

r 8 {M 9 + 23M 8 r + 233M 7 r 2 + 1365M 6 r 3 + 5098M 5 r 4 + 
12592MV + 20576M 3 r 6 + 21456M 2 r 7 + 12960Mr 8 + 3456r 9 } 
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For t — 1 the numerator becomes 

tt 2 M{768/(/ + l)M n + (12032/(7 + 1) + 1536)M 10 r + (76160/(/ + 1) + 24064)M 9 r 2 + 
(240384/(/ + 1) + 150272)M 8 r 3 + (350256/(7 + 1) + 142976)M 7 r 4 + 
(37968/(/ + 1) - 717248)M 6 r 5 - (469392/(/ + 1) + 436992)M 5 r 6 - 
(265040/(/ + 1) - 1506816)M 4 r 7 + (310048/(/ + 1) - 993792)M 3 r 8 + 
(124992/(/ + 1) - 884736)M 2 r 9 - (141696/(/ + 1) + 829440)Mr 10 + 27648r n } 

whereas the denominator reads 

/3 2 r 5 {16M 13 + 336M 12 r + 3016M n r 2 + 14928M 10 r 3 + 

43297M 9 r 4 + 69255M 8 r 5 + 37937M 7 r 6 - 52347M 6 r 7 - 81046M 5 r 8 + 

3504M 4 r 9 + 49376M 3 r 10 + 720MV 1 - 14688Mr 12 + 3456r 13 } 

For the last case, t — 2, the numerator is simply 

vr 4 M 3 (1024r- 128M) 

while the denominator, a bit more complicated, reads 

(3 4 r(r - 2M) 4 (2r + M) 2 (4r + M) 
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